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ABSTRACT 

Minimum  expected-cost  solutions  to  the  single-period 
inventory  problem  (the  Newsboy  problem)  have  been  well  known 
for  many  years,  for  risk  cases  where  the  distribution  of 
demand  is  known.  Also  well  known  are  minimax  cost  and 
Laplace  solutions  for  the  uncertainty  case  where  only  the 
range  of  demand  is  known.  This  study  explores  the  use  of 
order-statistic-based  quantile  estimators  as  decision  pro- 
cedures while  data  gathers  during  the  early  decision  periods. 
Simulation,  using  a variety  of  demand  distributions  and 
unit-cost  values,  provides  recommendations  on  which  period 
to  leave  the  minimax  procedure  for  a suggested  data-based 
decision  rule  which  will  be  asymtotic  to  the  optimal  risk 
procedure . 
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I.  INTRODUCTION 
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The  classic  single  period  inventory  problem  is  called  by 
various  well-known  names  such  as  "Newsboy  problem,"  "News- 
paperboy  problem,"  and  so  on.  [6].  In  this  thesis  it  will 
be  called  the  Newsboy  problem. 

Solutions  to  the  Newsboy  problem  are  well  known  under 
both  "risk"  and  "uncertainty"  conditions.  [7].  The  purpose 
of  this  thesis  is  to  study  the  Newsboy  problem  under  initial 
conditions  of  "Uncertainty,"  (i.e.,  when  the  possible  future 
for  demand  is  not  yet  known)  and  then,  for  the  subsequent 
periods  as  data  gathers,  to  suggest,  structure,  and  evaluate 
decision  procedures  during  this  transition  phase  to  decision 
under  "Risk." 

The  inventory  model  studied  in  this  thesis  is  a time 
independent  (stationary)  and  linear-cost  type.  The  study 
will  consider  the  problem  of  minimizing  the  expected  cost. 
Since  with  the  Newsboy  problem  both  the  expected  profit 
maximization  and  expected  cost  minimization  cases  have 
identical  solutions,  the  study  should  be  applicable  to  both 
interpretations  of  the  problem. 

Chapter  II  reviews  the  Newsboy  problem  for  decisions 
under  risk  and  decisions  under  uncertainty,  and  then  suggests 
possible  alternative  procedures  for  "from  uncertainty  to 
risk"  conditions.  Among  these.  Order  Statistics  (a  non- 
parametric  approach)  appear  promising  and  are  discussed  in 
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detail  in  Chapter  III.  Chapter  IV  evaluates  the  alternative 
nonparametric  transition  phase  procedures  of  the  Newsboy 
problem  and  a$  a result  proposes  a transition  phase  pro- 
cedure. In  Chapter  V conclusions  and  recommendations  for 
further  work  are  given. 


II . SOLUTION  OF  THE  SINGLE -PERIOD  INVENTORY  (NEWSBOY) 
PROBLEM  UNDER  RISK  AND  UNCERTAINTY  CONDITIONS 


In  this  section  the  Newsboy  problem  is  examined  under 
risk  and  uncertainty,  and  optimal  decision  rules  for  these 
conditions  are  given.  These  two  cases  form  the  bounds  of 
the  transition  phase  and  permit  us  to  characterize  the 
transition  phase  decision  problem  we  intend  to  study. 

Here,  when  it  is  said,  "under  risk,"  it  is  meant  that 
the  probability  distribution  of  demand  is  known.  When  the 
distribution  of  demand  is  not  known,  and  data  in  hand  is 
not  sufficient  to  estimate  the  possible  futures  of  demand, 
then  the  problem  is  called  "under  uncertainty."  [7], 


A.  SOLUTION  OF  THE  NEWSBOY  PROBLEM  UNDER  RISK 

Let  fD(x)  be  the  probability  density  function  of  the 
demand  for  continuous  populations,  and  p^(x)  be  the  proba- 
bility mass  function  for  discrete  populations.  If  CSQ  is 
the  stock-out  or  shortage  cost  per  unit  of  unsatisfied 
demand,  and  Cc  the  carrying  or  surplus  cost  per  unit  of 
excess  supply,  then  the  cost  for  a particular  period  of 
Newsboy  problem  will  be: 


Cost  C(S)  = 


|Cc(S-X) , if  0 < X < S,  and 


c cso CX-S) , if  X > s, 


where  X is  the  realization  of  demand  for  that  period,  and 
S is  the  inventory  level  in  hand.  Here  S represents  the 
decision  variable  for  this  problem. 
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The  expected  cost  E[C(S)]  which  results  from  keeping 
an  inventory  level  S in  hand,  will  be: 


E[C(S)]  = Z Cc(S-X)pn(X)  + Z C (X-S)  p (X) 
X=0  U X=S+lso  u 


for  discrete  demand  distributions  and. 


E[C(S)]  = Cc  /(S-x)  fD(x)dx  + CsQ  / (x-S)  fD(x)dx,  (2) 
0 S 

for  continuous  demand  distributions . 

Now,  the  expected  cost  for  a period  is  to  be  minimized 
and  we  wish  to  find  an  optimal  value  of  S,  say  SQ,  that  will 
yield  the  minimum  expected  cost.  In  other  words,  (1)  and 
(2)  are  to  be  minimized  with  respect  to  S. 

In  the  discrete  demand  distribution  case,  two  necessary 
conditions  for  a minimum  at  SQ  are: 


E[C(Sq  + 1)]  - E[C(Sq)]  > 0,  and 


E[C(S0-1)]  - E[C(S0)]  > 0. 

These  conditions  yield  the  well  known  result  [7] 


pD(So‘i:)  - Cc  + Cso  -pD(So)’ 


as  conditions  for  a minimum  expected-cost  solution  at  SQ. 

Similarly,  for  continuous  demand  distributions,  to 
minimize  (2),  the  first  derivative  of  E[C(S)]  is  taken  with 
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respect  to  S,  and  after  it  is  set  equal  to  zero,  yields 


VSo>  " 


so 


C + C 
c so 


(4) 


where  Fjj(Sq)  is  the  cumulative  distribution  function  (CDF) 
of  demand  at  SQ . Thus,  the  SQ  values  which  satisfy  equations 
(3)  and  (4)  are  the  optimum  inventory  levels  to  keep  on  hand 
in  order  to  minimize  the  expected  cost. 

This  solution  (under  risk)  is  quite  idealistic,  since  it 
assumes  that  the  distribution  of  demand  is  known.  In  fact, 
in  reality  we  are  seldom  certain  about  the  distribution.  in 
most  cases,  neither  the  distribution  of  demand,  nor  an  ade- 
quate amount  of  data  to  provide  a good  basis  for  estimation 
of  the  distribution  is  initially  available.  The  situation 
when : 

1.  Distribution  of  demand  is  not  known,  and 

2.  Enough  data  from  the  system  to  estimate  the  distri- 
bution of  demand  is  not  available,  but 

3.  We  do  have  some  idea  about  the  range  of  demand  D, 
say  0 < D < Dmax,  is  one  we  shall  call  "uncertainty." 

B.  SOLUTION  OF  NEWSBOY  PROBLEM  UNDER  UNCERTAINTY 

It  has  been  shown  that  whether  minimax  cost,  minimax 
regret,  or  the  Laplace  approach,  assuming  a uniform  distri- 
bution of  demand  is  chosen  as  a principle  of  choice,  the 
optimal  decision  under  uncertainty  is  the  same  [6],  [7], 
namely, 
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So  = 


C D 
so  max 

C + C 
c so 


, for  continuous  demand. 


So  < 


Cso(Dmax  + 1) 
Cc  + Cso 


< S0  + 


for  a discrete  demand  distribution.  The  results  of  these 
uncertainty  procedures  are  not  affected  by  newly  generated 
demand  data  unless  it  causes  one  to  revise  the  range  of  demand. 

C.  TRANSITION  PHASE  OF  THE  NEWSBOY  PROBLEM 

Now  the  two  bounds  or  limit  solutions  (uncertainty  and 
risk  conditions)  for  the  transition  phase  of  the  Newsboy 
problem  are  reviewed.  Implicitly,  to  proceed  in  accordance 
with  these  current  solutions,  one  begins  under  Uncertainty 
and  uses  the  Uncertainty  optimum  inventory  level  SQ  (based 
on  the  range  of  demand) , until  demand  data  is  adequate  or 
complete  to  estimate  the  distribution  of  demand  and  change 
to  the  optimum  inventory  level  SQ  for  the  Risk  conditions. 

As  mentioned  before,  the  purpose  of  this  study  is  to  find  a 
good  transition  phase  procedure  for  the  intermediate  periods 
when  some  data  has  accumulated  and  is  available.  A transition 
phase  procedure  may  be  characterized  by  two  decisions: 

1.  When  to  begin  the  transition  phase  (i.e.,  at  what 
period  should  we  leave  Uncertainty  procedures  to  begin  to 
use  an  estimator  of  S ) . 

2.  What  estimator  of  SQ  to  use. 

After  explaining  the  nature  of  the  transition  phase  in 
the  next  chapter,  three  candidate  estimators  for  the  optimum 
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inventory  level  SQ  to  use  during  transition  phase  will  be 
proposed. 

As  we  remember,  the  optimal  choice  of  S under  Risk  was 
based  on  a quantile.  Let  us  call, 


P 


so 


C + C 
c so 
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and  then  from  equation  (4)  we  have, 


s0  = xp  = fd' Cp> 


which  provides  a minimum  expected-cost  solution  under  Risk. 
Thus,  the  problem  is  to  locate  the  p quantile  of  the 
Cumulative  Distribution  Function,  F^(X) . As  more  and  more 
data  is  collected,  one  might  be  able  to  estimate, 


A 

s 

G 


~ - I 

F A(p)  , 


as  a function  of  the  data.  Some  alternative  approaches  to 
estimate  this  quantile  were  considered  to  be: 

1.  To  consider  similar  systems  for  which  we  have  demand 
data,  and  then  to  apply  a Bayesian  approach. 

2.  To  decide  about  the  form  of  the  distribution  of 
demand  and  carry  out  a maximum  likelihood  procedure  to  find 
the  parameters. 

3.  To  apply  a nonparametr ic  approach. 

For  the  first  two  cases,  if  one  is  able  to  decide  about 
the  form  of  distribution  such  as  lognormal,  gamma,  compound- 
poisson  process,  beta,  etc.,  then  the  p1"*1  quantile  of  demand 
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III.  A NONPARAMETRIC  APPROACH  FOR  THE  NEWSBOY 
PROBLEM,  TRANSITION  PHASE 

In  this  section  the  possibilities  for  applying  a non- 
parametric  or  distribution- free  approach  to  the  problem  of 
the  transition  phase  of  the  Newsboy  problem  will  be  explored. 
After  explaining  the  reasoning  for  a nonparametric  approach, 
three  candidate  estimators  of  the  optimum  inventory  level 
SQ  (or  X ) will  be  proposed. 

A.  WHY  A NONPARAMETRIC  APPROACH? 

As  mentioned  at  the  end  of  the  preceding  section,  if  the 
decision  maker  is  able  to  predict  the  form  of  the  demand 
distribution,  then  a parametric  approach  can  be  more  efficient 
and  have  a higher  convergence  rate  to  the  true  value  of  the 
estimator  *-han  a nonparametric  equivalent.  Hadley  and 
Whitin  [5],  give  various  considerations  relating  to  choice 
of  a demand  distribution  and  about  inventory  demand  predic- 
tions, in  detail.  For  a parametric  approach,  such  as  would 
result  from  fitting  a theoretical  distribution  to  data,  they 
suggest  to  use  poisson  distribution  when  the  lead  time  is 
relatively  low,  or  normal  distribution  when  it  is  not  rela- 
tively low,  etc. 

Since  it  is  quite  difficult  to  specify  the  form  of  the 
demand  distribution,  a nonparametric  approach  was  considered 
to  be  of  particular  interest  for  this  problem.  Also,  Hadley 
and  Whitin  point  out  that,  for  single-period  inventory  models. 


M 
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like  the  Newsboy  problem,  the  empirical  distribution  function 
(nonparametric)  is  reasonably  applicable,  unless  we  are  forced 
to  operate  the  model  at  the  tail  of  the  distribution.  [5]. 

Performance  of  the  model  decision  rules  was  assumed  to 
be  critical  during,  say,  the  first  50  periods  when  data  is 
inadequate.  The  decision  maker  needs  to  operate  the  system 
well  during  that  critical  time.  Later  on  he  will  have  a 
sufficient  amount  of  data  to  treat  the  problem  as  one  under 
Risk  and  simply  use  the  optimal  decision  rules  for  Risk.  For 
these  early  periods  (transition  phase)  a nonparametric  approach 
was  used. 

B.  THE  EFFECT  OF  THE  ESTIMATORS  OF  OPTIMUM 

INVENTORY  LEVEL  ON  EXPECTED  COST 

Before  we  suggest  some  estimators  of  optimal  inventory 
level,  Xp,  it  may  be  useful  to  see  how  those  estimators 
affect  the  total  expected-  cost  in  a period.  In  the  work 
that  follows,  it  is  assumed  that  the  demand  process  is  time 
independent  (stationary) , which  means  that  seasonal  effects 
and  trends  are  removed.  Also,  we  assume  that  for  the  previous 
period  we  know  what  the  demand  was. 

Before  we  continue  on  with  the  mathematical  development, 
it  is  useful  to  expand  on  the  conditions  for  an  optimal 
expected-cost  solution.  When  the  demand  distribution  is 
assumed  known,  we  remember  from  equation  (2)  that  the  expected 
cost  for  the  case  of  continuous  demand  is: 

S 00 

E[C(S)]  - Cc  / (S-x)  fD(x)dx-CS0  / (x  - S)  fD(x)dx. 


To  find  an  optimum  inventory  level  SQ,  taking 


E'  [CCS)  ] = d E|Cg(S)J 


S = S 


(5) 


and  setting  it  equal  to  zero  (necessary  conditions  for 
minimization) , we  obtain: 

- Cso  * Cso  FD<So>  * Cc  VSc>  - 0 • °r 


W ■ r ■ *°c  ■ p • and  so  ■ Ccp). 

c so 


Similarly, 

E~[C(S)] 


d2  E TCCSD 1 


S 


Cso  *WCc  W 

^Cso  + Cc^  W > 0 ’ 


(6) 


which  forms  a sufficient  condition  for  minimization.  Here, 

t"  H 

SQ  is  the  value  that  corresponds  to  the  p quantile  of  the 
population,  say  X . 

We  are  interested  in  the  case  where  the  population  is 
not  known,  but  there  exists  a sample  of  size  n of  demand 
data,  acquired  during  the  last  n periods. 

Let  Xp  X^,....,  Xn  be  realization  of  demand  during  the 
last  n periods.  Although  we  don't  know  the  value  of  Xp  yet, 
we  do  have  a chance  to  estimate  this  value  as  a function  of 
available  data,  say: 

Xp  * g(Xx,X2 Xn)  * Fp1  (p)  = So  (7) 
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To  express  the  expected  cost  as  the  function  of  Xp,  let 

5 = S - S = X - X . 
o o p p 

Equation  (2)  then  becomes: 

s + 5 

E[C(S0)]  = Cc  f°  (S 0*«-x)  fD(x)dx*Cso 


* Cso  's  ‘ So  ■ 3)  fD(x)dx 
°o 


If  equation  (8)  is  expanded  using  the  Taylor  series  around 


SQ,  it  gives 


E[C(Sq)]  = E [C  C S o ) ] + E'[C(Sq)]  E[5]  + E~[C(S0)]  ^ 


Since  E"[C(Sq)]  = 0 by  equation  (5),  and  having  E"[C(Sq)] 
from  equation  (6) , then  for  three  terms  the  Taylor  expansion 


Ef  ■521 


E[CCS0)]  . E[C(S0)]  * tCso*Cc)  fD(S0)  ^ 


Equation  (9)  shows  that: 

1.  Under  Uncertainty  conditions  the  expected  cost  for 
a period  is  linear  in  the  mean  square  error  of  the  pt^1 
quantile  estimators,  i.e.,  in 

E[62]  - E [ (Xp  - Xp)2] 

2.  If  one  is  able  to  find  a consistent  estimator 

A 

Xp(n),  of  Xp,  then  when  n gets  large  the  mean  square  error 
will  converge  to  zero  and  E[C(Sq)]  will  converge  to  E[C(S  )], 

A 

providing  that  the  moments  of  Xp(n)  exist  and  behave  properly. 
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3.  The  optimal  performance  for  an  estimator  here  is  to 

yield  the  expected  total  cost  under  Risk,  E[C(S  )]. 

Now,  we  wish  to  find  a good  consistent  point  estimator 

Xp  for  Xp.  Also,  it  is  not  required,  but  it  is  desirable  to 

obtain  one  that  is  an  unbiased  estimator  of  X . If  so, 

P 

A 2 A 

E[(Xp  - Xp)  ] will  be  equal  to  the  variance  of  estimator  Xp. 
Then,  by  equation  (9),  if  it  is  possible  to  improve  this 
estimator  reducing  its  variance,  it  will  directly  reflect 
and  will  decrease  the  expected  cost  accrued  by  using  that 
estimator . 

To  obtain  candidate  nonparametr ic  estimators  of  Xp, 

order  statistics  will  be  used.  In  the  following  section 

we  shall  develop  three  different  extimators  of  Xp , using 

an  order  statistics  approach.  Additional  details  about  order 

statistics  and  empirical  distribution  functions  related  to 

quantile  estimation  are  given  in  Appendix  A. 

1 . The  First  Candidate  Estimator  of 
Optimum  Inventory  Level:  X^ 

Let  Xn,  < X,1A  < ...  < X,  be  a set  of  n reali- 
(1)  (2)  (n) 

zations  from  an  unknown  distribution  in  order  statistics 

a l_  ^ L 

form,  where  X^  is  the  r largest  value  (r  order 
statistic).  Let  r be  defined  in  the  following  way: 

np  , when  np  is  integer 

r = 

[n  p + 1]  , otherwise , 

where  p is  the  quantile,  n is  the  number  of  periods  so  far, 

and  [X]  denotes  the  integer  part  of  X.  For  large  n,  X^ 

is  a consistant  and  unbiased  point  estimator  of  Xp  (Appendix  A). 
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In  Figure  1,  the  order  statistic  X^  is  represented 
graphically  with  the  distribution  function  FD(X)  and  the 
empirical  distribution  function  Gp(X) • 

Know  results  of  this  estimator  include  the  following. 
For  a uniform  (0,1)  distribution. 


E [X 


r 

n+1 


(10) 


and 

cov[Xf  ,,Xf  ] = cov[Uf  vUf  J = r(n-S-*-) , for  r < s 

irj  '•SJ  irj  lsJ  (n+1) z (n+2) 

(ID 

Similarly,  for  any  other  distribution. 


E[X(r)J 


F " 1 f -r-  -1 

rD  Ln+lJ 


(12) 


and 


cov[X 


y r (n-s  + 1)  . 1 

(r) ’ tS)  ' ^(Fd'1^)) 


provided  r < s. 

Later  these  results  will  be  used  to  compare  the 
efficiency  of  the  estimators  of  Xp. 


(13) 


2 . The  Second  Candidate  Estimator 

of  Optimum  Inventory  Level:  X(r) 

As  a second  candidate  estimator  of  Xp , X^  is 

offered  to  smooth  the  empirical  distribution  function  by  a 

linear  interpolation.  It  might  be  considered  that  the  area 

between  GD(X)  and  the  actual  distribution  function  FQ(X) 

may  be  a measure  of  deviation.  So,  here  interpolation  was 

done  at  the  midrange,  yielding  a new  empirical  distribution 
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function  R,  which  would  possibly  have  a smaller  area  between 
R and  FD(X),  than  the  area  between  GQ(X)  and  Fq(X) . A possibl 
realization  for  X^  is  shown  in  Figure  2.  Since  the  area 
between  those  distribution  functions  is  subject  to  change 
because  of  randomness,  the  results  for  that  case  would  be 
evaluated  numerically  rather  than  analytically. 

3 . The  Third  Candidate  Estimator 
of  Optimum  Inventory  Level: 

As  a third  candidate  estimator  of  X , a pooled  order 

P 

statistic,  combining  other  order  statistics,  is  proposed. 

For  such  an  approach  to  be  successful  by  equation  (9) , it 
has  to  reduce  the  mean  square  error  to  a lower  value  than 
the  rt^1  order  statistic  does,  if  it  is  possible  to  make  this 
pooled  estimator  an  unbiased  estimator  of  X^.  Then  the 
problem  is  one  of  comparing  the  variance  of  this  estimator 
with  the  variance  of  the  rt^1  order  statistic. 

To  examine  this  subject,  it  is  reasonable  to  begin 
with  the  uniform  (0,1)  distribution,  since  the  expectations, 
variances,  etc.  of  its  order  statistics  are  simple  and 
exact,  and  then  to  discuss  the  assumptions  that  are  necessary 
before  we  generalize  it  for  other  distributions, 
a.  ^(r)  for  t*ie  Uniform  Distribution 

Let  X , X,  X,,..  be  the  order  statistics 

(r-1)’  (r)  (r+1) 

from  the  uniform  (0,1)  distribution.  In  addition  to  the  r^ 
order  statistic,  a pooled  estimator  in  the  form  of  a linear 
combination  of  the  (r-1),  (r) , (r+1)  order  statistics  will 
be  considered.  Here  our  third  candidate  estimator  is: 


FIGURE  2.  Second  Candidate  Estimator 


X (r)  ' Xp  = a X (r-1)  + 3 X(r)  + Y X(r+1)  » 
where  a + 8 + y = 1. 

As  mentioned  earlier,  if  this  estimator  can  be  an  unbiased 

estimator  of  X , then  a comparison  of  the  variances  of  the 
1*  Vi 

r order  statistic  (First  candidate  estimator),  and  this 
pooled  estimator  (Third  candidate)  will  indicate  which  one 
has  better  performance. 

To  maintain  unbiasedness,  we  must  have: 


E[Xfr)]  - E[X  ,] 


(r)J  ’ 


E[X(r)]  ■ a E[X(r_i;)]  * 6 E[X(r)]  * Y E[X(r+1)] 


If  demand  has  a uniform  distribution,  then 


EtX(r)l  " * + 6 + Y thP  ' 

This  quantity  is  equal  to  E[X^r^]  = only  when  a = y. 

One  way  to  decide  on  the  values  of  a,  6 and  y is  to  choose 
values  so  that  the  variance  is  minimized.  So,  we  would  like 
to  minimize 

VarC?(r))  - vartc  X^.^S  X(r)  * Y X(r+1))  , 
subject  to. 


a + 3 + y = 1 


and  a = y 


By  the  constraints  on  a,  3 and  y,  these  coefficients  can  be 
written  in  terms  of  a only,  yielding 


Var(X(r))  - Var{a(X(r.1}  * X(rn))  * (l-2«)  X(r)> 
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or, 

VarCX(r))  = a2(Var(XCr.1))  + Var(X(r+1))) 

+ 2 Cov(XCr-l)’X(r+l)) 

+ (2cc-4ct2)  (CovCX^.^.X^)  + Cov(X(r3,X(r  + 1))) 
+ Cl  - 4a  + 4a2)Var(X(r))  .. 

Taking  the  derivative  of  VarfX^.^)  with  respect  to  a and  setting 
it  equal  to  zero  gives, 

a = 2 Var(Xtr1)  - CovCXfr-11  .Xfrp  - Cov(Xfr1  .Xfr4,1-|) 

Var(X(r-i,)  *Var(XCrtl))  *2Cov(X(r.1),XCrn))  - 

4[Cov(Xr.q>X(r))  *Cov(X(r),X(rtl))  Var(X(r))  ] 

(14) 

Having  available  the  values  of  variance  and  covariances  by 
equation  (11),  equation  (14)  becomes 

2r(n-r+l)  - (r~l)(n-r+l)  - r(n-r) 

a 

(r-1) (n-r+2)  + 4r(n-r+l)  + (r+1) (n-r)  - 4(r-l) (n-r+1)  - 

4r{n-r)  + 2(r-l)(n-r) 

the  n,r  terms  in  the  numerator  and  in  the  denominator  cancel, 
resulting  in  a = 1/2.  So,  to  minimize  the  variance  of  X ^ , 
the  coefficients  have  to  be  a = y = 1/2,  and  3=0. 

This  says  that  when  a population  is  uniform,  a third 
reasonable  estimator  of  Xp  is: 

?(r)  ’ ^V-l)  * x(r.l)>  - ‘15> 
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yielding 

Var(f(r) 


2r(n-r+lj  - (n+1) 

2(n+l)2  (n+2) 


If  the  ratio  cf  the  two  candidate  estimators  is  taken, 
2r (n-r+1) 


Var(X(r)) 


Var 


CX(0) 


2r(n-r+l)  - (n+1) 


always.  This  implies  that  for  a uniform  distribution, 

Var  (X^)  > Var  (X^r^),  and  is  the  better  estimator  of 


Xp.  In  Figure  3,  for  n - 20,  these  variances  are  sketched. 


When  r is  close  to  bounds  (either  close  to  1 or  n) , the 
difference  between  those  two  variances  become  larger. 


r) 


FIGURE  3.  Variances  of  X.  . and  X.  . for  Uniform(0,l) 
Distribution , when  n is  twenty  . 
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i 


These  results  apply,  of  course,  only  to  uniform 


demand  distribution.  In  the  following  paragraphs  the  gener- 
alizations to  any  demand  distribution  will  be  discussed. 

b.  Generalization  of  X^  to  any  Demand  Distribution 

For  distributions  other  than  uniform,  it  is 
difficult  to  claim  that  , as  given  by  equation  (15),  is 

exactly  unbiased.  Also,  the  variance  and  covariances,  as 
equation  (13),  will  be  functions  of  the  density  fQ(X).  Even 
though  these  two  problems  exist,  some  reasonable  approxima- 
tions can  be  discussed.  If  we  have  both  a large  n and  fairly 
flat  densities,  since  (r-1),  (r) , (r+1)  are  the  closest 
neighbors  among  the  other  order  statistic  values,  the  densi- 
ties corresponding  to  these  points  would  tend  to  equal  each 
other  and  during  the  calculation  of  a,  as  in  equation  (14), 
would  essentially  cancel. 

If  we  have  flat  densities  (where  the  three  points 
are  nearly  linear),  or  if  we  have  a large  number  of  observa- 
tions (when  the  points  are  close  together  accordingly), 

X(r)  would  be  approximately  unbiased,  too.  Thus,  it  may  be 
reasonable  in  some  cases  to  use  X^  for  other  distributions 
besides  the  uniform. 

Incidentally,  we  would  not  like  to  try  another 
pooled  estimator  of  X such  as 

*(r)  a2  X (r=2) +al  X(r-1)  + B X(r)+Yl  X(r+1) + Y2  X(r+2) 
For  this  estimator,  the  assumptions  which  were  made  to  permit 


■ 
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X^r-j  to  apply  the  distributions  other  than  the  uniform  may 
no. longer  be  appropriate.  For  example,  in  equation  (13), 

fD*-FD  1 (n+T-^  and  fD^FD  1 would  probably  have  quite 

different  values,  because  these  are  the  densities  of  the  two 

farthest  order  statistics,  Xf  and  Xr  in  X,  . . 

I r-  l ) (^r+z ) t r ) 

Otherwise,  if  these  densities  were  known,  we  would  not  use 
a nonparametric  approach. 

So  far,  solutions  of  the  Newsboy  problem  under 
"Uncertainty"  and  under  "Risk"  have  been  reviewed,  and 
possible  procedures  for  a transition  phase  between  them 
discussed.  Also  in  this  chapter,  three  candidate  estimators 
of  the  optimal  inventory  level  Sq  (or  X ) were  developed  for 
use  as  alternative  transition  phase  procedures.  In  the  next 
chapter,  the  performance  of  these  alternative  transition 
procedures  will  be  evaluated  by  a computer  simulation  of 
Newsboy  problem  decision  making.  Finally,  from  the  simula- 
tion we  will  obtain  information  about  what  time  to  leave  the 
minimax  (Uncertainty)  rule,  for  each  transition  phase  pro- 
cedure . 


IV.  EVALUATION  OF  THE  CANDIDATE  ESTIMATORS 
OF  OPTIMUM  INVENTORY  LEVEL,  FOR  THE 
TRANSITION  PHASE  OF  THE  NEWSBOY  PROBLEM 

In  Chapter  III,  the  rationale  for  using  a nonparametr ic 

approach  for  the  transition  phase  of  the  Newsboy  problem  was 

discussed.  It  was  shown  that  under  Uncertainty  the  expected 

cost  for  a period  was  approximately  linear  in  the  mean  square 

error  of  the  optimum  inventory  level  SQ  (or  Xp)  corresponding 

to  the  pt^1  quantile,  and  three  candidate  estimators  of  Xp  or 

S were  introduced, 
o 

In  this  chapter,  the  performance  of  these  estimators 
will  be  tested  for  several  demand  distributions  by  simulation 
After  explaining  our  reasons  for  using  simulation,  we  will 
explain  the  structure  of  the  simulation  program  used.  In 
part  C,  the  simulation  parameters  used  in  the  program  are 
described,  together  with  the  reasons  for  choosing  these 
parameters.  Part  D illustrates  the  results  of  simulation. 
These  results  will  be  used  later  to  compare  the  performance 
of  the  estimators  of  Xp  and  the  transition  phase  procedures, 
and  to  see  whether  there  were  any  apparent  adverse  effects 
of  approximations  or  assumptions. 

A.  RATIONALE  FOR  SIMULATION 

The  purpose  of  the  study  was  to  find  a transition  phase 
procedure  (i.e.,  which  estimator  of  Xp  or  SQ  to  use  and  when) 
during  the  transition  from  Uncertainty  to  conditions  under 
Risk  for  the  Newsboy  problem.  As  explained  before,  the 
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performance  of  an  alternative  transition  phase  procedure 
depends  directly  on  the  performance  of  the  candidate  esti- 
mator of  Xp  used,  in  that  there  is  a one-to-one  correspondence. 

A primary  reason  for  simulation,  of  course,  was  to  study 
the  question  of  how  long  one  should  continue  using  the  mini- 
max procedure  (for  Uncertainty)  before  changing  to  a transi- 
tion phase  using  a quantile  estimate.  Simulation  over  several 
demand  distributions  should  provide  insight  regarding  when 
(with  respect  to  cost)  this  change  in  the  Newsboy  problem 
should  occur.  Additionally,  we  wish  to  further  examine  the 
question  of  which  of  the  candidate  estimators  to  use. 

However,  it  was  hoped  that  the  estimators  of  X or  S 

p o 

would  come  close  to  having  desired  order  statistics  proper- 
ties, being  exact  for  the  case  of  the  uniform  distribution 
(equations  (9),  (10),  (11)  and  (12)).  For  example,  for  the 

rt^1  order  statistic  X,  , it  has  been  shown  that 

(r) 

Fp  (E  [X  £r  -j  ] ) = for  a uniform  distribution,  but 

FD  (e  [x  (r)  ^ ^ for  convex,  and  FQ  (E  [X  (r -j  ] ) >_  for 

concave  cumulative  distribution  functions  of  demand.  [4]. 
Unfortunately,  there  is  not  an  analytical  measure  to  show 
how  much  the  slack  in  these  inequalities  would  change 

\ ; 

depending  upon  the  amount  of  convexity  or  concavity  of  the 
demand  distributions. 

Also,  for  the  third  candidate  estimator  of  X , some 
additional  assumptions  were  made  that  might  not  hold  for 
every  distribution. 


i 
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In  summary,  our  reasons  for  simulation  are: 

1.  To  find  the  proper  time  (period)  to  begin  the 
transition  phase  (to  use  one  of  the  candidate  estimators), 
leaving  the  "Uncertainty"  case  and  the  minimax  decision  rule, 

2.  To  provide  information  to  supplement  incomplete 
analytical  knowledge  of  order  statistics  regarding  popula- 
tion quantiles, 

3.  To  evaluate  the  impact  of  the  approximations  made 

in  developing  the  third  candidate  estimator,  , over 

various  demand  distributions. 


B.  DESIGN  FOR  SIMULATION 

The  simulation  program  was  written  in  the  FORTRAN  IV 
language.  In  every  simulation  run,  inventory  stocking 
decisions  are  made  successively  using  one  particular  demand 
distribution  and  one  particular  value  of  the  quantile 

C0 

(P  = — ) , say,  demands  generated  from  a normal  parent 

for  quantile  0.3  ...  etc.  The  demand  distributions  and 
various  quantile  values  used  in  the  simulation  are  given  in 
the  next  section.  Each  run  consisted  of  50  successive 
decision  periods  starting  with  no  demand  information  other 
than  the  range  estimate,  and  replicated  60  times. 

As  output  values  of  the  program,  costs  for  each  period 
were  averaged  over  the  number  of  replications.  These  average 
costs  per  period  were  calculated  for  five  different  cases: 

1.  The  Uncertainty  case,  with  the  minimax  decision  rule 
being  used,  denoted  by  M. 

2.  The  ideal  case  (Risk),  with  the  quantile  value  known, 
denoted  by  ID. 
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3.  The  transition  phase  case  using  estimator  Xr  , denoted 

by  R.  ('rJ 

4.  The_transition  phase  case  using  estimator  X , , denoted 

by  R.  1 

The  minimax  (for  Uncertainty)  and  ideal  (for  Risk)  cases 
would  give  constant  inventory  levels,  i.e.,  the  same  amounts 
all  the  way  through  the  simulation  run.  For  Uncertainty,  the 
range  of  demand  was  approximately  known  and  would  not  change 
unless  it  was  revised,  while  for  Risk  conditions  it  was 
assumed  that  the  distribution  of  demand,  and  thus  the  quan- 
tile, was  known. 

Transition  phase  optimum  inventory  level  values  generated 
by  X^  , and  X ^ ^ were  subject  to  change  from  period  to 

period  until  finally  they  converge  to  the  value  associated 
with  the  Ideal  case.  For  the  first  two  periods,  ^(r) ’ 

would  not  lead  to  an  inventory  level,  since  they  need 
at  least  two  data  points. 

The  program  was  run  for  50  periods,  and  replicated  60 
times.  So,  period  costs  for  each  transition  phase  procedure 
were  averaged  over  60  replications.  Program  steps  can  be 
ordered  in  the  following  way: 

1.  Define  the  demand  distribution  yielding  the  quantile 
p,  give  the  stockout  cost  CSQ,  the  carrying  cost  CQ,  and  the 
range  of  demand. 

2.  Compute  and  set  minimax  and  ideal  case  optimum 
inventory  values. 

3.  Generate  a random  demand  value  from  the  specified 
demand  distribution. 
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4.  Find  the  single-period  cost  for  each  case  (minimax, 
X(r)’  ^(r)’  *(r)  an<*  Ideal)  by  comparing  their  inventory 
levels  and  the  demand  generated  for  this  period.  Record 
these  costs.  Find  the  next  period  inventory  levels  revising 

X(r)’  ^(r)  ’ ^(r)  anc*  t0  t*ie  next  period. 

5.  Repeat  steps  3 and  4 for  50  periods. 

6.  Again,  beginning  from  the  first  period,  repeat 
step  5 for  60  times  (replications) , and  at  the  end  of  60 
replications  find  the  average  costs,  and  the  sample  variances 
of  costs  over  60  replications,  for  each  period. 

7.  Tabulate  the  results. 

The  results  obtained  from  this  simulation  will  be  examined 
in  part  D. 

C.  SIMULATION  PARAMETERS 

In  this  simulation  a variety  of  possible  forms  of  distri- 
butions, including  continuous,  discrete,  symmetric,  skewed 
(long-tailed) , etc.,  were  chosen.  Demand  distributions  and 
their  parameters  used  in  this  simulation  were  the  following: 

Distribution  Parameters 


Uniform 

Poisson 

Chi-square 

Lognormal 


the  range  over  [0,30] 

with  mean  10  units 

with  10  degrees  of  freedom 

generated  by  a normal  parent 
having  a mean  of  3 and  a variance 
of  1 

with  mean  35  and  variance  100 
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Normal 


For  D in  the  Minimax  procedure,  y + 3t  is  used  for 
max 

''SO 

tailed  distributions.  Quantiles  p = ^ ^ — , tested  were 
0.10,  0.30,  0.50,  0.70,  0.90,  0.95  for  each  distribution. 

These  were  determined  so  that  CsQ  + Cc  was  constant  (in  this 
simulation  CSQ  + Cc  = 4.0  was  used).  So,  for  each  {quantile, 
distribution}  combination  the  differences  of  the  expected 
period  costs  would  come  from  the  performance  of  the  estimators 
of  Xp,  as  was  shown  earlier  by  equation  (9). 

In  the  next  section  we  shall  present  and  discuss  some 
of  the  results  from  the  simulation. 

D.  SIMULATION  RESULTS 

In  this  part,  simulation  results  will  be  evaluated  and, 
for  the  Newsboy  problem,  a transition  phase  procedure  will 
be  proposed.  In  presenting  these  results  we  will  identify 
the  quantile  value  as  p.  Tables  I and  II  give  overall  mean 
costs  and  standard  deviations  of  cost.  Then,  Table  III  will 
show  at  what  period  the  transition  procedure  should  begin 
for  a particular  distribution,  and  quantile.  Table  IV  will 
give  relative  efficiency  of  the  best  transition  procedure 
to  the  Ideal  case  (decision  under  Risk) . 

On  the  next  page,  Table  I gives  the  overall  mean  values 
of  the  cost  per  period  (over  50  periods  and  60  replications) 
for  each  transition  case,  using  CsQ  + Cc  = 4.00.  As  can  be 
seen  from  Table  I,  there  does  not  seem  a significant  differ- 
ence using  R,  R,  or  R,  even  though  R looks  quite  consistent 
at  the  tail  ( p=  0.95).  One  reason  may  be  the  fact  that  the 
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TABLE  I 


Average  Costs  Over  50  Periods  and  60  Replications 
having  CQ  + CsQ  = 4.0 


M=Minimax,  R=X^r^,  R=X ^ r ^ > R=Xj-rj,  ID=Ideal  Cases 


p 

Uniform 

Poisson 

Chi-Square 

Lognormal 

Normal 

M 

5.49 

3.23 

3.01 

67.00 

11.66 

R 

6.35 

2.38 

2.91 

14.50 

8.09 

0.1 

R 

6.31 

2.38 

2.88 

14.50 

8.09 

R 

6.30 

2.38 

2.87 

14.48 

8.07 

ID 

5.49 

2.04 

2.56 

11.44 

6.93 

M 

12.68 

4.60 

5.92 

255.96 

20.04 

R 

13.58 

4.55 

5.98 

36.34 

14.25 

0.3 

R 

13.55 

4.52 

5.96 

36.23 

14.18 

If 

13.46 

4.56 

5.94 

36.38 

14.16 

ID 

12.68 

4.26 

5.62 

30.15 

13.17 

M 

15.08 

5.97 

11.08 

339.86 

15.71 

R 

16.  03 

5.35 

7.47 

51.22 

15.80 

0.5 

R 

15.97 

5.33 

7.45 

51.37 

15.78 

If 

15.91 

5.34 

7.45 

51.35 

15.73 

ID 

15.08 

4.98 

6.98 

43.49 

14.80 

M 

12.75 

7.33 

11.81 

298.67 

14.97 

R 

13.75 

4.74 

7.21 

57.24 

14.24 

0.7 

If 

13.66 

4.74 

7.21 

57.23 

14.19 

If 

13.60 

4.72 

7.20 

57.66 

14.19 

ID 

12.75 

4.44 

6.67 

48  72 

13.24 

M 

5.39 

4.37 

6.26 

132.42 

9.33 

R 

6.31 

2.68 

4.40 

42.36 

7.76 

0.9 

R 

6.30 

2.67 

4.38 

42.38 

7.76 

R 

6.30 

2.66 

4.36 

42.58 

7.73 

ID 

5.39 

2.33 

3.71 

35.91 

6.70 

M 

2.83 

2.38 

3.34 

70.42 

5.27 

R 

3.77 

1.75 

3.00 

31.80 

5.04 

0.95 

R 

3.76 

1.75 

2.99 

31.74 

5.02 

t 

3.76 

1.75 

3.00 

32.44 

5.04 

ID 

2.83 

1.41 

2.27 

25.47 

3.91 
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difference  between  variance  of  and  variance  of  , 

was  not  much  (Figure  3).  But  if  one  wishes  to  see  the  dif- 
ference no  matter  how  much  it  is,  the  following  normalization 
by  ideal  costs  may  be  considered: 

1.  For  a distribution  from  Table  I,  for  every  quantile 
R.-ID 

£ find  — — , and  sum  this  value  over  all  quantiles  g.  for 
every  R^,  where  R^  is  R,  R or 

2.  Repeat  step  1 for  every  distribution. 

3.  Pick  the  R^'s  which  have  the  smallest  sum  among 
R,  R,  R for  every  distribution  as  their  best  procedure. 

The  results  of  this  normalization  were  as  follows: 


Uniform 

Poisson 

Chi-square 

Lognormal 

Normal 

R 

0.871 

0.767 

0.859 

1.253 

0.839 

R 

0.845 

0.  752 

0.831 

1.251 

0.824 

f 

0.828 

0.755 

0.821 

1.295 

0.816 

The  best  transition  procedures  were  to  use  R for  uniform, 
E for  poisson,  R for  chi-square,  R for  lognormal,  and  E for 
normal  distributions.  To  see  the  effects  of  the  very  right 
tail  of  a distribution  on  transition  procedures,  p = 0.99 
was  used  for  normal  and  chi-square  distributions.  At  this 
quantile  the  minimax  costs  were  below  R,  R,  and  R costs. 

This  says  minimax  can  be  used  at  the  very  tail  of  a distri- 
bution along  the  transition  phase,  providing  very  accurate 
prediction  of  the  range  of  demand  (say  . Otherwise,  if 

we  remember  from  Table  I for  p = 0.95,  R,  R and  1 were 
better  than  minimax  approach. 


From  Table  II  (overall  standard  deviations  over  50  periods 
and  60  replications),  it  is  seen  again  that  there  is  not  much 
difference  of  standard  deviations,  even  though  R has  a better 
trend  (smaller  sample  variance  of  cost)  among  R,  JT  and  R. 

From  these  two  tables  it  can  be  said  that  there  is  not 
much  difference  of  the  costs  yielded  using  R,  R and  R,  but 
if  an  order  is  necessary  it  should  be  the  one  using  R for  p 
less  than  0.95,  and  R at  the  tail  (i.e.,  when  p is  greater 
than  0.95) . 

During  the  simulation,  as  the  number  of  periods  increase 
the  estimators  converge  to  Xp  or  the  true  SQ,  while  the 
minimax  and  ideal  case  inventory  levels  were  constant. 

Since  the  values  of  X^,  X^  , X^  usually  change  from 
period  to  period  because  they  are  updated  by  the  most  recent 
demand  value,  they  will  generally  be  associated  with  a larger 
variance  of  cost  than  the  minimix  and  ideal  cases.  This  is 
reflected  in  Table  II.  Sometimes  the  minimax  procedure  had 
a smaller  sample  variance  of  cost  (while  possibly  yielding 
the  largest  average  cost)  than  the  ideal  case.  This  may  be 
explained  as  follows: 

The  ideal  solution  we  are  using  relates  to  minimizing 
expected  cost,  not  the  variance  of  cost.  These  two  minimi- 
zation points  are  generally  different,  depending  on  the 
distribution  function.  If  the  minimax  solution  is  different 
from  the  ideal  (expected  cost)  inventory  level  and  closer  to 
the  minimum  variance  solution,  it  would  yield  a smaller 
variance  of  cost.  Further  derivations  about  the  variance 
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TABLE  II 


Sample  Standard  Deviations  over  50  Periods 
and  60  Replications 

M=Minimax,  R=Xj-r-j,  R=Xj-r ^ , R^X^  , ID=Ideal  Cases 


p 

Uniform 

Poisson 

Chi-Square 

Lognormal 

Normal 

M 

3.15 

1.26 

1.79 

43.46 

3.  77 

R 

6.36 

2.69 

3.27 

27.46 

9.89 

0.1 

R 

6.26 

2.65 

3.21 

27.47 

9.70 

R 

6.26 

2.59 

3.18 

27.47 

9.53 

ID 

3.15 

1.62 

2.26 

15.90 

6.62 

M 

7.46 

3.18 

4.61 

76.10 

10.33 

R 

9.00 

3.60 

4.93 

57.64 

12.17 

0.3 

R 

8.84 

3.53 

4.89 

57.69 

11.97 

K 

8.85 

3.51 

4.86 

57.71 

11.99 

ID 

7.46 

3.16 

4.42 

45.93 

10.95 

M 

8.82 

4.30 

6.35 

66.95 

12.34 

R 

10.31 

4.14 

6.37 

84. 29 

12.  73 

0.  5 

R 

10.23 

4.12 

6.33 

83.91 

12.68 

R 

10.23 

4.14 

6.44 

84.23 

12.75 

ID 

8.82 

3.84 

5.80 

72.72 

12.09 

M 

7.30 

3.47 

4.87 

46.79 

9.62 

R 

9.08 

3.95 

6.66 

99.68 

11.65 

0.7 

I 

9.01 

4.04 

6.83 

99.68 

11.65 

R 

8.87 

4.09 

6.76 

99.23 

11.67 

ID 

7.30 

3.55 

5.66 

92.87 

10.42 

M 

3.16 

1.25 

1.87 

32.45 

3.90 

R 

6.59 

3.28 

5.72 

99.99 

9.10 

0.9 

R 

6.55 

3.33 

5.76 

99.23 

9.09 

R 

6.  55 

3.29 

5.69 

98.05 

9.01 

ID 

3.19 

2.57 

4.10 

94.42 

6.34 

M 

1.66 

0.69 

1.12 

31.66 

2.07 

R 

6.50 

2.91 

5.36 

94.64 

8.50 

0.95 

R 

6.50 

2.95 

5.34 

92.85 

8.47 

R 

6.50 

2.94 

5.23 

92.74 

8.46 

ID 

1.66 

1.49 

3.15 

84.97 

4.11 
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minimization  solutions  of  the  Newsboy  problem  are  given  in 
Appendix  B. 


The  simulation  results  about  the  period  to  begin  a 

transition  phase  are  given  in  Table  III.  Shown  are  the 

first  periods  in  which  the  average  period  cost  (over  60 

replications)  using  the  estimator  was  below  that  obtained 

using  the  minimax  rule.  In  the  simulation  it  was  observed 

that  once  costs  using  an  estimator  fell  below  the  minimax 

costs,  they  stayed  below  for  the  remaining  time  periods. 

Simulation  results  showed  that,  on  the  average,  at  the  fifth 

period  we  were  supposed  to  begin  a transition  procedure, 

ceasing  to  use  the  minimax  procedure.  In  other  words,  for 

the  first  five  periods  the  large  variability  of  candidate 

estimators  of  SQ  (or  X ) caused  a larger  average  cost  than 

the  minimax  procedure  yielded.  As  a very  special  case,  if 

the  distribution  of  demand  is  uniform,  minimax  and  Ideal 

case  solutions  would  match  each  other  and  thus  the  uniform 

distribution  is  omitted  from  Table  III.  Also  by  coincidence, 
Cso 

if  p = ~ — — t* — is  at  or  close  to  an  intersection  of  the 
Lc  Lso 

implied  uniform  CDF  of  minimax  procedure  and  the  actual 
distribution  function  of  demand,  then  the  minimax  procedure 
may  yield  the  Ideal  case  solution.  Otherwise,  minimax  in- 
ventory level  has  always  a higher  cost  than  the  Ideal  case. 

In  Table  III  the  large  numbers  are  those  specific  values 
which  relate  to  a quantile  P that  is  close  to  minimax-- 
actual  CDF  intersection  point.  So  it  takes  time  for  R,  R 
and  R to  yield  costs  less  than  the  minimax.  No  matter  how 


TABLE  III 


Periods  to  Begin  Transition  Phase 
Leaving  the  Minimax  Procedure 


juantile  (p)  Poisson  Chi-Square  Lognormal  Normal 


0.95 


* Ideal-Minimax  intercept  points  (or,  close  to  these  points) 


At 


long  it  will  take  to  fall  below  the  rainimax  cost  (here, 
nearly  at  period  ten)  from  simulation  outputs  it  was  observed 
that  after  the  fifth  or  sixth  period  the  cost  difference 
between  the  minimax  and  the  transition  procedures  was  not 
large.  This  also  can  be  predicted  from  the  pattern  values 
of  Table  III.  The  numbers  in  Table  III  are  relevant  for  all 
candidate  estimators,  i.e.,  R,  £ and  £ all  fell  below  the 
minimax  cost  simultaneously  at  the  same  period. 

As  a measure  of  efficiency  for  a transition  procedure, 
the  ratio  of  average  cost  induced  by  this  procedure  to  Ideal 
average  cost  is  used.  Table  IV  gives  these  ratios  for  the 
average  costs  over  second  and  fifth  10  periods  for  R (since 
R seems  the  promising  transition  procedure  among  the  other 
ones)  . 

We  will  examine  the  data  of  Table  IV  for  the  second  10 
period,  assuming  that  the  minimax  has  been  left  and  the 
transition  phase  has  begun.  We  also  believe  that  convergence 
toward  the  true  value  of  Ideal  inventory  level  SQ  for  the 
third,  fourth,  ...  10  periods  will  not  be  worse  than  the 
convergence  accomplished  during  the  second  10  period. 

If  the  transition -procedure  (here,  R)  was  operated  at 
the  tail  (i.e.,  p >_  0.9),  then  the  efficiency  ratio  for  the 
second  10  periods  was  on  the  average  1.16  and  more  (worse). 

As  the  quantile  gets  close  to  1.0,  the  minimax  procedure  was 
yielding  much  better  results  than  a transition  procedure 
(assuming  the  range  of  demand  is  predicted  quite  well) . 
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TABLE  IV 


The  Ratio  of  the  Second  and  Fifth  Ten-Period 
(11  to  20,  41  to  40)  Average  Costs  £o  the 
Ideal  Case  (Under  Risk)  Costs,  Using  X.  . 

Transition  Procedure  ^ 


Quantile 

(P) 

Uniform 

Poisson 

Chi-Square 

Lognormal 

Normal 

0.1 

S 

F 

1.063 

1.000 

1.065 

1.023 

1.058 

1.046 

1.033 

1.002 

1.108 

1.008 

S 

1.065 

1.059 

1.050 

1.040 

1.066 

0.3 

F 

1.007 

1.041 

1.003 

1.004 

1.008 

S 

1.056 

1.069 

1.039 

1.041 

1.077 

0.5 

F 

1.013 

1.021 

1.013 

1.004 

1.007 

S 

1.073 

1.069 

1.061 

1.057 

1.054 

0.7 

F 

1.000 

1.001 

1.006 

1.006 

1.004 

S 

1.102 

1.113 

1.128 

1.190 

1.096 

0.9 

F 

1.008 

1.016 

1.052 

1.026 

1.027 

S 

1.170 

1.193 

1.211 

1.240 

1.183 

0.95 

F 

1.002 

1.040 

1.054 

1.062 

1.055 

Notes : 

S stands  for  "SECOND,”  F stands  for  "FIFTH." 

When  p is  less  than  0.9  overall  ratio  is  1.060  for  second 
and  1.010  for  fifth  10  periods. 
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For  the  rest  of  the  quantiles  where  p is  less  than  0.9, 

this  ratio  for  the  transition  procedure  was  almost  constant 

1.060.  To  test  the  effects  of  quantiles  and  the  effects  of 

the  form  of  demand  distributions  on  the  efficiency  ratio 

(1.060),  a two-way  Analysis  of  Variance  could  be  applied  to 

the  data  of  Table  IV  for  quantiles  less  than  0.9. 

Different  distributions  were  significant  at  a level  0.03, 

and  the  different  quant'iles  were  significant  at  a level  of 

0.01.  This  result  was  explained  that  for  C + C =4.0 

so  c 

cost  units,  whatever  the  distribution  of  demand  was  and  for 
any  quantile  less  than  0.9,  the  ratio  of  the  cost  induced  by 
transition  phase  procedure  (X^  , here)  to  the  Ideal  cost 
(expected  cost  under  Risk)  would  be  practically  constant  at 
1.060  for  second  10  period,  and  would  not  exceed  this  amount 
for  the  further  periods  (for  example,  at  the  fifth  10  period 
this  ratio  was  1.01). 

On  the  next  part,  as  the  summary  of  this  chapter,  and 
using  the  results  obtained  in  it,  a Single-Period  inventory 
(Newsboy)  problem  transition  phase  procedure  is  suggested. 

E.  PROPOSED  NEWSBOY  TRANSITION  PROCEDURE  ALGORITHM 

(1)  Approximate  the  bounds  of  Demand  (Range). 

(2)  Set  p = c C+s°  . 

(3)  If  p > 0.90,  go  to  step  6,  otherwise  go  to  step  4. 

(4)  Use  minimax  for  the  first  5 or  6 periods. 

(5)  For  6 <_  n <_  50,  use  R (unless  a reason  for  R,  E)  , 
and  go  to  step  10. 

(6)  If  p <_  0.95,  go  to  step  8,  otherwise  go  to  step  7. 
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(7)  If  step  1 is  accurage,  may  use  minimax  for  all  periods 
and  go  to  step  10,  but  if  it  is  not  accurate,  go  to 
step  8. 

(8)  Use  minimax  for  the  first  5 or  6 periods,  then  go  to 
step  9. 

(9)  For  6 <_  n <_  50,  use  R and  go  to  10. 

(10)  Analyze  data. 

*May  fit  a probability  distribution  to  the  data  and  begin 
Decision  Under  Risk  case. 


*May  use  R,  R or  R for  n > 50  too.  Since  they  converge 
X quite  fast,  the  risk  of  using  them  after  50  periods  hope- 
fully will  be  nearly  equal  to  the  risk  of  Ideal  case. 
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r.  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 

For  the  Newsboy  problem,  a nonparametric  decision  pro- 
cedure that  allows  transition  from  Uncertainty  to  decision 
under  Risk  has  been  developed.  Three  estimators  of  the 
demand  corresponding  to  the  p1"*1  quantile  were  obtained, 

Cs0 

where  p = ~ — —7? — . Testing  these  estimators  over  several 
'■'C  + '-so 

distributions  of  demand  with  Newsboy  cost  criteria,  it  was 
found  that  excluding  cases  such  as  drastic  skewness,  pooling 
the  (r-l)st  and  (r+l)st  order  statistics  for  R provided 
essentially  the  best  estimator,  where  R and  R also  worked 
well . 

It  was  also  found  that  for  the  distributions  tested  and 

Q 

excluding  (r--  - s,°r  ■)  values  very  close  to  0.0  or  1.0,  that 
'-so  uc 

the  switching  Newsboy  decision  making  from  minimax  solutions 
to  quantile  estimator  could  be  made,  possibly  at  the  sixth 
period. 

For  further  work,  some  parametric  studies  such  as  log- 
normal distribution,  or  beta  and  gamma  families,  depending 
upon  whether  the  population  of  demand  is  finite  or  not,  may 
be  tried.  Then,  given  that  the  decision  maker  does  not  know 
what  the  distribution  of  demand  is,  which  procedure,  nonpara- 
metric or  parametric,  gives  better  response  to  Uncertainty 
may  be  compared.  As  another  recommendation,  some  work  can 
be  afforded  to  improve  nonparametric  approaches  further. 
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APPENDIX  A 


ESTIMATION  FOR  POPULATION  QUANTILES 
BY  . ORDER  STATISTICS 


In  this  Appendix  order  statistics,  especially  those 
concerning  population  quantiles,  will  be  briefly  reviewed. 
For  a further  study  of  order  statistics,  references  [1], 
[2]  and  [3]  are  cited. 

The  Joint  Distribution  of  n Order  Statistics 


Order  statistics  theory  is  mainly  based  on  continuous 
distribution  functions.  In  reference  [2]  for  discrete 
parents  the  joint  distribution  of  order  statistics  is  derived 
as  well. 

Let  X^,X2>--.,Xn  be  a random  sample  from  a probability 
distribution  function  F^.  If  X^^,i=l,2,...,n  is  the  i1"*1 
largest  value  in  this  sample,  then  X^  < X^j  < • • • < ^(n) 
are  the  order  statistics  of  a random  sample  X^  ,X2 » . • • ,Xn . 

1*  U 

Here,  X,  . is  called  the  r order  statistic.  Since  Fv  is 
(r)  X 

a distribution  function,  by  the  probability  integral  trans- 
formation U = Fx(x)  will  have  uniform  distribution  over  [0,1], 


l . e . , 


FyCu)  = 


u,  for  0 < u < 1 
0,  otherwise. 


Letting  Fx(x)  = u with  Fx  (X)  £ u if  and  only  if  X < x, 
the  following  equations  hold: 

FyCu)  = P(U<u)  = P{Fx(X)<u}  * P [X<_FX_1  (u)  ] = Fx[F*1(u)]=u 
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Since  F^(X)  is  a monotone  increasing  function  and  by 
equation  (16)  there  is  a one-to-one  correspondence  between 


Fx(x(i))  < Fx^(2)  < •••<Fx(X(n)^  and  X(l)  < X (2) 


Let  fx  be  the  density  of  F^,  so  that 


f y y y (x, ,x7 , . . . .xl  is  the  joint  density  of 

A1,x2,...,xn  i i n 


the  sample  ,X2 , . . . ,Xn<  If  each  X^,  i=l,2,...,n 

is  independent  and  identically  distributed,  then 


n 


■ Y Y Y ^ 1 > ? 

a1,a2* • * * 9 An  L 


(xn  ,X,f.  . . ,Xn)  = TT  f Y (x  • ) . 


i = l 


X^i- 


Since  there  exists  n!  permutation  of  those  sample  values, 


n 


;x(1),x(2) x^OV^-'-'V  = • (17) 


is  the  joint  density  of  n order  statistics. 

From  the  previous  references  cited: 

fX(r-)  = FI)  !-{iTTTT  [FX(yr)ir  rfx(yr)" 

(18) 

When  F^  is  a uniform  distribution  over  (0,1),  this  equation 
turns  to  be. 


n!  r-1  n-r 

£X(r/n^  = (r-1)  ! (n-r)  ! ^ [1'ul 


(19) 


and  the  following  sequence  of  equations  are  given. 


E‘X(r)kl  * (T-linn-D!  ' ^ 


and 
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t 


y (r+k-1) (r+k- 2) . . . . (r+1) r 
E^X(r)  ^ ~ (n+k) (n+k- 1) . . . . (n+2) (n+1) 

I 

where  k is  a positive  integer.  For  k = 1 : 

EfXCr)l  ‘ HTT  ' 

By  equations  (20)  and  (21) 

Var[X(r)]  ■ E[X(r)2]  - [E(X(r))l2  - 

and , 

E [X  r . - X,  .]  = — rlg  + U , 

1 (r)  (s)J  (n+l)(n+2) 


(20) 


(21) 


(22) 


(25) 


From  equations  (6)  and  (8), 

C°v[X(r),x(s)l  * E(X(r)X(s)]  - E[X(r)lE[X(s)l,  and 


Cov  [X(r),X(s)] 


r (n-s+1) 
(n+1)2  (n+2) 


(24) 


Now,  we  wish  a transition  from  the  uniform  distribution 
to  a general  distribution  F^.  Using  the  relations  in 
equation  (16) , 


U(r)  = Fx^X(r)^  » or  X(r)  = FX  (U(r))  » 

where  U,  . is  the  rth  order  statistic  from  Uniform  (0,1). 

(r) 

Then,  E[X(r)]  - E^'1  (U(r))]  * 1 [E(U(r))]  = Fx" 1 (jjJj)  . 

(25) 


Furthermore,  the' references  given  at  the  beginning  of  this 
appendix  show  that,  for  continuous  distributions  having 
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= [npj  + 1] , and  0 < px  < p2  . . . < pA  < 1 , 

the  joint  distribution  of  X(ni) >X(n2) ’ ' * * X(n£)  asymtotically 
tends  to  an  ^-dimensional  normal  distribution  with  means 
p^,  p2,...,  p^  , and  covariances. 


Mi-PiJ 

cov[X(n  ,,X(  ,]  = n f fx  ) f (x  ) , provided  i < k 

1 K Pi  Pk 


(26) 


So  the  variance  of  X , . will  be, 

(n^) 


Var[X,„  ,]  - 


PjCl-Pi) 


fni)'  n{f(Xp  )}2 


Having  n^  = r,  it  can  be  proven  that, 


U = 


X(ni)  'XPi 


^ <x(rrV 


/Var (X(n. ) } /P~TrrpT7 


fx(Xp.)  % N(0,1)  (27) 


Equation  (27)  is  quite  important  for  this  study.  It 
says  that,  whatever  the  distribution  function  is,  the  rt}l 
order  statistic,  X(r)>  approximately  normally  distributed 
around  (distribution  free  property).  This  also  means 
that  asymtotically  X^  is  an  unbiased  estimator  of  Xp . 

Let  Gq(X)  be  the  empirical  distribution  function  of 


order  statistics  X(]j»X(2) 

9 • • • 

>X(n)  having 

1 

r 0 , 

| 

if 

X < x(l) 

GDCX)  - < 

| ” ’ 

if 

X(J)  i x * x 

1 

1 > 

if 

X * X(n) 

■1) 
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In  Reference  [3]  it  is  said,  by  the  Glivenko-Centelli 
theorem  that,  GQ(X)  converges  uniformly  to  FX(X)  and  nGQ(X) 
is  binomial  (n,Fx(X)),  i.e., 

p[nGD(X)-j]  = p[GD(X)  = j]  = (")[Fx(X)]J[l-Fx(X)]n‘j, 

(28) 

Point  and  Interval  Estimation  of  Population  Quantiles 

For  point  estimation  of  a population  quantile,  we  assume 
£ h 

that  Xp  is  the  p quantile  of  the  distribution  function  Fx, 
and  P (X  < Xp)  = Fx(Xp)  = p , or  Fx" 1 (p)  = Xp. 

a.  l_ 

Then,  by  equation  (27)  the  r order  statistic  is  a consistent 

and  unbiased  estimator  of  the  p**1  quantile,  where 

” np  , if  np  is  integer 
r = ■ C 

„ [np  + lj  , otherwise  , and  p = ^ 

Lso  Lc 

For  interval  estimation  of  population  quantiles,  we  are 
looking  for  r and  s such  that, 

P^X(r)  < Xp  < Xj-s-j)  = y , where  y is  percent  confidence. 

Now, 

P<X(r)  < Xp  < X(s)>  = P(X(r)  < V ' P(X(s)<  V«  or 
by  equation  (16) 

P(X(r)  < Xp  < X(s)>  ■ < P)  - PCFX<X(s)>  < P)- 

Since  p = FY(X  ),  by  equation  (28) 
a p 

P(X(r)  <Xp<X(s))  = Jr(  i )Pi(1"Pf1-  S ( Pi(1-Pf"i 

s- 1 

- £ (?  )p1(l-pf 1 = Y (29) 

r 
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To  summarize  these  derivations: 

1.  For  large  n the  r1"*1  order  statistic  is  a consistent, 

^ h 

unbiased  point  estimator  of  the  p quantile  for  any  distri- 
bution. 

t*  Vi 

2.  Confidence  intervals  for  the  p quantile  as  given 
by  equation  (29),  distributed  by  the  Binomial  distribution 
(n,p)  . 

3.  As  n gets  larger,  X^  converges  to  and  also  the 
empirical  distribution  function  GQ(X)  converges  to  actual 
cumulative  distribution  function  FX(X) . 
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APPENDIX  B 


EXPECTED  COST  AND  VARIANCE  MINIMIZATION  POINTS 
OF  THE  NEWSBOY  PROBLEM 

In  this  Appendix  we  will  attempt  to  explain  the  differ- 
ence between  minimization  of  the  expected  cost,  and  mini- 
mization of  the  variance  of  the  cost,  for  the  Newsboy 
problem. 

Decision  makers'  utility  functions  are  involved  with 
both  expected  cost  (or  return) , and  the  variance  of  this 
cost.  If  a decision  maker  is  risk  prone,  the  variance  of 
the  cost  may  not  affect  him,  and  he  just  tries  to  minimize 
expected  cost.  On  the  other  hand,  if  he  is  risk  adverse, 
then  he  may  give  up  some  of  the  cost  minimization  value 
for  a smaller  variance  (a  trade-off  between  expected  cost 
and  variance  of  cost).  As  an  example  [5],  let 

2 

U(x)  = x - bx  , be  a utility  function.  Then, 

E [U (x) ] = E [x]  - bE  [x2]  = E [x]  - b(Var(x)  + [E(x)]2). 

If  we  try  to  optimize  this  value,  we  not  only  deal  with 
expected  value  of  x,  but  also  variance  of  x. 

On  the  following  pages,  the  variance  minimization  points 
of  some  distributions  related  with  the  Newsboy  problem  are 


derived . 
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When  this  is  applied,  to  a uniform  distribution  over 


c°-  W-  i-6-’  fD(x)  = 


D ’ 0 i x i Dmax  and  FD(S> 

max 


then,  the  left  hand  side  (LHS)  of  equation  gives, 

LHS  = S2(C  2 - C 2 ) — + C 2 S - C 2 -5^*  , 

c so  ZD  so  so  Z 

max 


Similarly,  right  hand  side  (RHS)  of  it  yields, 

Rhs  • (%^)s-cS0)(#--  (cc  + c ) -c  s + 

max  max 

After  we  equate  LHS  - RHS  and  opened  them,  it  was  found  that, 

C D 

S = s° max  , (32) 

0 C + C 

c so 

is  the  point  that  minimizes  the  variance  of  the  cost  for 
the  uniform  distribution.  This  point  is  exactly  the  same 
point  that  minimizes  expected  cost. 

Let  us  apply  the  same  procedure  to  a bell-shape,  say, 

Beta  (a  = 2,  B = 2).  Then 

£D(X)  = r (aVr'(g)  Cl-*)6'1  = 6x(l-x)  for  0 < x < 1, 


Fn(x)  = f 6u(l-u)du  = 3x2  - 2x3,  and 
u 0 


1 - Fd(x)  = 1 - 3x2  + 2x3,  I u fD(u)du  = y(l-4x3  + 3x4) , 


/ uf.  (u)du  = y (4x3  - 3x4  ) . 

0 U Z 

Putting  those  values  in  equation  (2),  gives  a 5 degree 
polynomial : 
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0 


■w 


re  + c )2s5-i(c  +c  )2s4  + 3(c  +c  )2  s3  + 

v c so'  2V  c so'  c so 


+ 4(C2  + C2  + 5C  C )S2  - (4C  2 + 5CC  )s4(C  C + 2C  2 ) = 
2V  c so  c so'  ^ so  c so'  2 ^ c so  so' 

(33) 

It  may  be  difficult  to  solve  this  equation,  but  it  is 
not  difficult  to  see  that  it  doesn't  have  the  same  solution 

Wlth  T^Tr — = 3 S 2 - 2S  3 , or  3S  2 - 2S  3 - p = 0, 

Cc+Cso  0 0 0 o 

which  is  the  minimum  expected  cost  solution. 

In  the  exponential  {A}  case, 

fD(x)  = Ae'Xx,  Fd(x)  = l-e'Xx,  l-FD(x)  = e‘Xx  and 


D' 


/ u f(u)du  = f[l-e'Xx(l+  x)], 
0 . A 


/ u f(u)du  = i[e'Xx(l+Ax.)] 
x A 


Again,  inserting  these  values  into  equation  (2)  gives, 

fd(S)  cca 


UTfe  » 

c so 


where  Fp(S)  = 1-e  X^. 


(34) 


If  we  compare  (34)  with  Fp(So)  = — 

P C S 0 

1 p 

or  SQ  = • yin  (^-  +-£ — ) , we  will  see  that  they  are  different 
c so 

points.  If  we  examine  equation  (34)  further,  it  can  be 
shown  that  S > SQ,  i.e.,  the  variance  minimization  point  S 
is  closer  to  the  tail  than  the  minimum  expected  cost  point  S 


-=  l-e"XSo 
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